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■ Abstract 

> ; 

On ! We study homogeneous and isotropic cosmologies in a Weyl space- 

£^ | time. We show that for homogeneous and isotropic spacetimes, the 

^NJ ■ field equations can be reduced to the Einstein equations with a two- 

fluid source. We write the equations as a two-dimensional dynamical 
| system and analyze the qualitative, asymptotic behavior of the mod- 

""^ ■ els. We examine the possibility that in certain theories the Weyl 

q-i 1-form may give rise to a late accelerated expansion of the Universe 

^ | and conclude that such behaviour is not met as a generic feature of 

the simplest cosmologies. 



1 Introduction 

Recent observations regarding the evolution of the Universe have led sev- 
eral authors to propose a variety of departures from the standard cosmology. 
These proposals can be roughly grouped into two categories. First, there 
exists a dark energy of unknown nature which is responsible for the accel- 
erating expansion of the Universe, (cf. PjJ |2| for comprehensive reviews and 
references). Alternatively general relativity requires a modification at cosmo- 
logical distance scales [Hill! (for a pedagogical review see jH]). Less explored 
is the idea that the geometry of spacetime is not the so far assumed Lorentz 
geometry (see for example [0]). Due to its simplicity, Weyl geometry is consid- 
ered as the most natural candidate for extending the Riemannian structure. 
In this geometry the covariant derivative of the metric tensor is not zero. For 
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example, in the spirit of the idea of creation of the Universe from quantum 
fluctuations of "nothing", Novello et al [7j considered perturbations of the 
geometry of the form 

5 (V P 57^) = (SQp) g^, 

where Q p is the gradient of a scalar field. Other motivations for reconsid- 
ering Weyl geometry include the hope to remove the cosmological big bang 
singularity. 

Weyl geometry can be incorporated a priori in a theory or a posteriori, for 
example as a consequence of the variational principle involved. This is the 
case of the application of the Palatini method to gravitational Lagrangians of 
the form L = f (R) [HIE], where it turns out that Q p = d p In /' (R). However, 
the Palatini device suffers from serious problems and leads to inconsistencies 
when applied to general Lagrangians for the construction of a gravity theory 
(see for example PHIE] for a thorough critique of the Palatini device). In jT2] 
it was shown that a consistent way to incorporate an arbitrary connection 
into the dynamics of a gravity theory is the so-called constrained variational 
principle. 

In this paper we study homogeneous and isotropic cosmologies in a Weyl 
framework. More precisely, we explore the field equations derived with the 
constrained variational principle from the Lagrangian L = R + L matt er, under 
the condition that the geometry be Weylian |12j . If the matter Lagrangian is 
chosen so that ordinary matter is described by a perfect fluid with equation 
of state p 2 = (72 — 1) P2, we show that for homogeneous and isotropic space- 
times, the field equations reduce to the Einstein equations with a two-fluid 
source. The first fluid with equation of state pi = p\ (i.e. stiff matter), stems 
from the Weyl vector field and the second fluid describes ordinary matter as 
mentioned above. 

The plan of the paper is the following. In Section 2 we express the field 
equations in terms of the Levi-Civita connection and show that homogeneous 
and isotropic models can be interpreted as two-fluid models. In Section 3, 
following the method of Coley and Wainwright in jlSj, we write the equations 
as a two-dimensional dynamical system and analyze the asymptotic behav- 
ior of the models. In Section 4 we consider the two-fluid model resulting 
from a modification of the Einstein-Hilbert Lagrangian, cf. f|17|) . and study 
the possibility of providing a mechanism of accelerating expansion. Readers 
unfamiliar with Weyl geometry may find in the Appendix an exposition of 
the techniques involved. In this Appendix we also derive the Bianchi identi- 
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ties and write the Einstein tensor and the Bianchi identities in terms of the 
Levi-Civita connection. 

2 Field equations 

We recall that a Weyl space is a manifold endowed with a metric g and a 
linear symmetric connection V which are interrelated via 

Vp9ap = -Qn9ap, (1) 

where the 1-form is customarily called Weyl covariant vector field (see the 
Appendix for a full explanation of the notation involved below). We denote 
by D the Levi-Civita connection of the metric g a p. 

In ^2] it was shown that application of the constrained variational prin- 
ciple to general Lagrangians of the form L — f (R) , in the context of Weyl 
geometry, yields the field equations obtained via the metric variation in Rie- 
mannian spaces with a source tensor depending on the Weyl vector field. 
The simplest theory that can be constructed with the constrained varia- 
tional principle is obtained from the Lagrangian L = R. The field equations 
are (see [T2] ) 

G iH = -V (M Q V) + Q^Qu + g»u (V a Q Q - Q a Q a ) =: M^. (2) 

If we express the tensor M^ u in terms of the quantities formed with the Levi- 
Civita connection D and take into account of (|A.29|) . the field equations 
become 




In the case of integrable Weyl geometry, i.e., when = d^, the source 
term is that of a massless scalar field. Taking the divergence of (J3j) and using 
the Bianchi identities (1A.31J) we conclude that 

D^Qn = 0. (4) 

In this paper we will be concerned with spatially homogeneous and isotropic 
spacetimes. Therefore we have to make the assumption that Q M is hyper- 
surface orthogonal. That means that is proportional to the unit timelike 
vector field n M which is orthogonal to the homogeneous hypersurfaces, 

Q»=:qn^ Q 2 = Q„Q» = -q 2 . 
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Formally the field equations (J3J) can be rewritten as 

o 

G^ u = (pi + Pi) n^n u + pig^, (5) 

with 

3 

Pi = Pi = 4<? 2 , (6) 

and we see that the equation of state of the q— fluid corresponds to stiff 
matter. For spatially homogeneous and isotropic models, the field equations 
((HJ) become the system of the equations: the Friedmann equation 1 

d\ 2 k 1 o 
a) + a*=l q > 

and the Raychaudhuri equation 

'a ( a\ 2 k 3 , 

2 -+ " + - = — A ( l ■ 
a \a J a z 4 



Equation (jlj becomes 
which implies that 



a 

q-3-q = 0, 

a 



C 



where C is a constant. Therefore the energy density and pressure of the 
q— fluid evolve as a" 6 . 

In the following we assume that ordinary matter is described by a perfect 
fluid with energy-momentum tensor, 

T^v = (p 2 + P2) u^u u + p 2 g^u, (7) 

where denotes the fluid velocity. To preserve the homogeneity and isotropy 
of the spacetime it is necessary that 



Q" = qu> t , Q 2 = Q^ = -q 
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1 We adopt the metric and curvature conventions of Here, a (t) is the scale factor, 
an overdot denotes differentiation with respect to time t, and units have been chosen so 
that c = 1 = 8ttG. 
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Therefore we are dealing with a two-fluid model with total energy density 
and pressure given by 



P = pl+P2, P = Pl+P2, 



(8) 



respectively, where 

Pi = Pi, P2 = (72 - 1)P2, 
i.e., 71 = 2 and 72 < 71. The field equations take the final form 



(9) 




(10) 



- = -^[4 Pl + (3 7 2 - 2)p 2 ] 



(11) 



Pi 




(12) 



P2 = -3 7 2P2-- 

a 



(13) 



It is well known (see ^1]) that the field equations for one fluid can be written 
as a two-dimensional dynamical system for the Hubble variable H := a/ a and 
the density parameter Q := p/3H 2 . A drawback of this analysis is that it does 
not give a complete description of the evolution for closed models. In fact, 
at the time of maximum expansion the Hubble parameter becomes zero and 
therefore, the time coordinate defined in jHj p. 58, cannot be used past 
the instant of maximum expansion. Instead, we deal with closed models by 
defining the compactified density parameter u, (see [T5] ) 



We see that u is bounded at the instant of maximum expansion (H = 0) and 
also as p — >• 0, in ever-expanding models. 

3 Phase plane analysis 

We now adopt the Coley and Wainwright formalism for a general model with 
two fluids with variable equations of state (cf. ^3])- Assuming that 71 > 72 
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(14) 



tan 2 



or 
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and 71 > 2/3, we define the transition variable x [ — lj 1] 

P2 + Pi 

which describes which fluid is dominant dynamically. The total density pa- 
rameter Q = Qi + Q 2 can be compactified as in (J 14)) and the evolution 
equation of the variable x is obtained by applying the conservation equation 
to pi and p 2 . Defining a new time variable r by 

d 3 (7x - 72) fp 1 d 



dt 2 V 3 cos a; dr 

and with the same kind of manipulations as in jT3|, one obtains the following 
dynamical system 

duj 1 

— = — - (0 — x) cos 2a; cos 
ar 2 

^ = (1- X 2 )sin„, (16) 



where the constant b is 

, 3 (71 + 72) - 4 



> -1. 



3 (71 - 72) 
In our case, we always have 71 = 2. 

The phase space of the two-dimensional system (|TK|) is the closed rectangle 

D=[-7r/2,7r/2] x [-1,1] 

in the a; — x plane (see Figure 1). 

By inspection we can see that the line segment {(lj, x) £ D : uj = is 
an invariant set of ([16)1 . It consists of three trajectories, the line segment 

{(w,x) £ D : uj = tt/4, -1< x < 1} 

and the equilibrium points (vr/4, —1) and (7r/4, +1) . Similarly we can specify 
the following invariant sets. 

u) = — 7r/2 contracting empty models Q = 0, H < 

X = — 1 one- fluid models ^2 = 

X = +1 one- fluid models Q± = 

u; = 7r/4 expanding flat models Q = 1, H > 

uj = — 7r/4 contracting flat models f2 = 1, < 

uj = tt/2 expanding empty models Q = 0, H > 
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X 

A 



F.7 +1 



k = -1 
< 1 



k = +1 



n > i 



o < 1 



■fii = 



-tt/2 



-tt/4 



tt/4 



tt/2 to 



n 
i 



= i 



■n 2 = o 



H < 



Fi + M+ 
H > 



Figure 1: The invariant sets and equilibrium points of (fTfiJl 



It is easy to verify that the equilibrium points lie at the intersection of these 
sets. The equilibrium points customarily denoted as F± 2 , F-f 2 , M + , M~ are 



F r,2 

M + 
M~ 



to = tt/4: expanding flat model Q = 1, k — 0, H > 

a; = — 7r/4 contracting flat model Q = 1, k — 0, H < 

uj = tc/2 expanding Milne model f2 = 0, = —1, H > 

= —tt/2 contracting Milne model Q = 0, = —1, H < 



and the subscripts indicate which fluid survives. Each of these eight equilib- 
rium points corresponds to an exact solution of the Einstein equations. 

We shall carry out in some detail the stability analysis of the equilib- 
rium points of (fTfij) in the case 72 = 1, corresponding to dust. It turns out 
that linearization of is sufficient to determine the global phase portrait 
of the system. In fact, the derivative matrix J (u,x) of the vector field 
of (pH)j) is non-singular and, J computed at each of the eight equilibrium 
points, has two real eigenvalues. Therefore, the Hartman-Grobman theo- 
rem applies in the case of It is easy to verify that J computed at 
all equilibrium points is diagonal. Therefore, we conclude in a straightfor- 
ward manner that (—tt/2, —1) , (— 7r/4, +1) , (vr/4, +1) , (tt/2, —1) are saddle 
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points, (— 7r/2, +1) , (7r/4, —1) are unstable nodes and (— 7r/4, —1) , (7r/2, +1) , 
are stable nodes. The phase portrait is shown in Figure 2. 




Figure 2: The phase portrait of (fTBj) with 72 = 1. 

Regions III and IV correspond to expanding models. The is a past 
attractor of all models with Q > 0, i.e., the evolution near the big bang is 
approximated by the flat FRW model where the Weyl fluid dominates. Open 
models expand indefinitely and the evolution is approximated by the Milne 
universe at late time. Flat models expand indefinitely and the evolution is 
approximated by the flat FRW universe at late time. In both cases the "real" 
second fluid dominates at late times while the q— fluid becomes insignificant. 
On the other hand, any initially expanding closed model in region III, how- 
ever close to F 2 + , eventually recollapses and the evolution is approximated 
by the flat FRW model where the Weyl fluid dominates. 

In the case 72 = corresponding to a positive cosmological constant, a 
new equilibrium point (0,2/3) appears, denoted by E (see Figure 3). It cor- 
responds to the Einstein static model. Consider for example a closed model 
in region III starting close to the model. Its trajectory passes close to the 
Einstein model, indicating a phase of halted expansion and asymptotically 
approaches the de Sitter model _F 2 + . Similarly, an open model in region IV 
starting close to F^ asymptotically approaches the de Sitter model. This 
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Figure 3: The phase portrait of (fTSJl with 72 = 0. 

attracting property of the de Sitter solution for all expanding models is not 
restricted only to isotropic cosmology. In fact, the cosmic no-hair conjec- 
ture states that all expanding universe models with a positive cosmological 
constant, asymptotically approach the de Sitter solution (see for example 

We therefore conclude that the Weyl fluid has significant contribution 
only near the cosmological singularities. In expanding models the "real" 
fluid always dominates at late times and therefore the contribution of the 
Weyl fluid to the total energy-momentum tensor is important only at early 
times. 



4 Extensions 

The field equations (JSJ) constitute the generalization of the Einstein equa- 
tions in a Weyl spacetime in the sense that they come from the Lagrangian 
L = R. There is however an alternative view, namely that the pair (Q,g) 
which defines the Weyl spacetime also enters into the gravitational theory 
and therefore, the field Q must be contained in the Lagrangian indepen- 
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dently from g. In the case of integrable Weyl geometry, i.e. when = <9 M 
where is a scalar field, the pair (0, g^ u ) constitute the set of fundamental 
geometrical variables. A simple Lagrangian involving this set is given by 

L = R + ^Q„ (17) 

where £ is a constant. Motivations for considering theory ()17|) see can be 
found in [7J [T7] (see also [TH] for a multidimensional approach and [TH] for 
an extension of (|17|) to include an exponential potential function of 0). By 
varying the action corresponding to (fT7j) with respect to both g^ v and one 
obtains 

G^u — -^2~~ \QnQv - ^Q 2 9tiuj > and D0 = 0. (18) 

Note that the second of (|18|) comes from the variational procedure while (J3J) 
is a consequence of the Bianchi identities. 

For isotropic cosmologies one can again interpret the source term in (|18|) 
as a perfect fluid with density and pressure given by 

Pl = Pl = \q 2 , A = ^^ (19) 

respectively. However, if we allow for £ to be a free parameter, equation (fT9"j) 
implies that the energy density and pressure of the q— fluid may take negative 
values. An important consequence is that open models (k = —1) in vacuum 
with A < avoid the initial singularity [7]. 

In the following we assume that A < and apply again the two-fluid 
analysis of Section 3. A sufficient negative contribution of the q— fluid to 
the total energy density and pressure, cf. (jSJ), may provoke an accelerating 
expansion of the universe. In fact, the p + 3p term in the Raychaudhuri 
equation 

a 1 . 

- = -- (p + 3p ), 

a o 

becomes negative provided that 4pi < (2 — 372) p2- Unfortunately, this can- 
not explain the observed acceleration of the Universe, as we shall see in a 
moment. 

The Friedmann constraint (fTU|) implies that the total density parameter 
Q = Qi + Q 2 is still non-negative for flat and closed models, but Q may be 
negative for open models. We conclude that Q cannot be compactified as in 
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(|14|) and the transition variable x defined by (J15)) is unbounded. Therefore 
we cannot apply the Coley and Wainwright formalism developed in Section 
4. 

Nevertheless, we can infer about the asymptotic behavior of the system 
by looking at the original field equations, (fTUjl - ffT^j) . From these equations 
we see that the state (a, a, pi, p 2 ) £ IR 4 of the system lies on the hypersurface 
defined by the constraint (jl(J|) and the remaining evolution equations can be 
written as a constrained four- dimensional dynamical system. By standard 
arguments one can show that for flat and open models the sign of H is 
invariant (see for example [20]), therefore an initially expanding universe 
remains ever expanding. Suppose that a (to) > at some time to- Then no 
solution of (fTTj) -(|T3 |) exists such that a (t) > for all t > t . In fact, (fT2|) and 
(fT3"j) can be solved to give 

* = ^P * = 5R5T. (2°) 

where C\ and C 2 are constants. Since 7 2 < 2, the "real" second fluid in 
an expanding universe eventually dominates. Therefore, the term 4pi + 
(372 — 2) p 2 in (fTTj) becomes positive at some time t\ > t and evidently 
a (t) < for all t > t±. We conclude that even if the universe expands initially 
with acceleration, eventually evolves according to the Friedmann cosmology. 

Note that for flat and open models there exist solutions without an initial 
singularity. To see this, consider for example 72 = 1 and substitute (|2Uj) into 
the Friedmann equation (jlOj) with k = — 1. Then (fTU|) can be interpreted as 
describing the motion of a particle with total energy E — 1/2 in the effective 
potential 

V» = -- + 4, A,B>0. (21) 
a a 4 

In Figure 4, we see that motion is impossible for values of the scale factor 
smaller than some minimum a m \ n and therefore, a (t) > a m ; n for all t. A 
similar argument shows that flat models, also avoid the initial singularity, a 
conclusion already obtained using phase portrait analysis by Oliveira et al 

To summarize, the extension to Weyl geometry in theory (fTTj) cannot ex- 
plain the observed acceleration of the Universe. Since the real fluid dominates 
at late times, the accelerated expansion due to the Weyl fluid is important 
only at early times. Nevertheless, it is possible that theory (fT7j) could pro- 
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Figure 4: The effective potential (|2*Tj). 

vide a geometric explanation of an inffationary phase present in the early 
universe. 
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Appendix: Weyl geometry 

The connection V of the Weyl space is determined by the symmetric functions 

-pa fa \ | f*ta. 

where {^ 7 } are the Levi-Civita connection coefficients and the tensor field 
C| 7 is given by 

C/3 7 = g ( S pQi + s "Qp ~ 9pyQ a ) ■ 

Given a metric and a connection satisfying ((H), viz. V ^g a p = —Q^gap, then 
for every differentiable function a, the metric and the 1-form defined by 

g' a p = vg«p, Q f ft = Qv- % h ° (A.22) 
12 



respectively, also satisfy (JTJ). Thus, g a p and are far from unique; rather 
g belongs to an equivalence class [g] of metrics and for each g E [g] , there 
exists a unique 1-form Q such that (0) is satisfied. A particular choice of a 
pair (Q, g) is called a gauge and (jA.22|) is a gauge transformation. 
The Riemann tensor is defined by 

R a p 1 s = <9 7 r /3<5 _ ^% + T^Tpg - T^ s Tp r (A. 23) 

In contrast to the familiar property in Riemannian geometry, in Weyl geom- 
etry is not antisymmetric in its first two indices; it satisfies 

2R(pv)af3 = g/iiyH al3 , (A. 24) 

where 

H a0 := d a Q p - dpQ a . (A.25) 

Note that H a p is the same for all derivative operators, for example H a p = 
V Q <5/3 — V^Qcj. When first introduced by Weyl, H a p was supposed to rep- 
resent the Faraday 2-form. The relation 2R^ a/3 = g^ v H a ^ can be obtained 
either directly from (jA.23|) . or more quickly by applying on the metric the 
anticommutativity property of V Q : 

(VqV/j — V pV a ) g^y = R f + R P u/3a9fj,p- 

The Ricci tensor is defined by 

Raj3 — R^ctpP 

and is not a symmetric tensor as we are accustomed to in Riemannian ge- 
ometry; its antisymmetric part is given by 

R[a/3] = H a p. 

The geometric meaning of (IA~24|) is that the measuring units change from 
point to point (see [2T] for a discussion). In fact, the length I 2 = g l _ lu V fl V 1 ' 
of a vector field V is not invariant when V is parallely translated in a small 
circuit ; instead, using the standard argument one may show that for a closed 
circuit enclosing an elementary area SS a/3 , the variation of the length is 

5 (I 2 ) = R 0a , )a pV>'V6S t # = l 2 H af3 6S a P. 
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As a consequence, there is an additional loss of synchronization for two ini- 
tially synchronized clocks following different paths from one point to another, 
due to the distinct variation of the units of measure along the two paths. This 
is the so-called second clock effect |2~2"] . 

When the Weyl vector field is the gradient of a scalar function, then 
the curl H a p vanishes identically and we have the so-called integrable Weyl 
geometry. In that case the spacetime is not a genuine Weyl space, but a 
conformally equivalent Riemann space. In fact, if = d^, where is a 
scalar field, it is easy to see that it can be gauged away by the conformal 
transformation g a p = (exp 0) g a p and therefore the original space is not a 
general Weyl space, but a Riemann space with an undetermined gauge 
However, most studies of gravity theories were developed in the framework 
of integrable Weyl geometry in which the second clock effect is eliminated 

m hhi na nni 

Bianchi identities 

For the convenience of the reader we present a few identities. We denote by 
D the Levi-Civita connection. 

1. V a Q a + Q 2 = V a Q a 

2. □ := g a/3 V a Vf3 = V Q V Q = V Q V a - Q a V a 

3. Q a V a = Q a V a 

4. = D^A" + C; Q A«, V^B U = D^B U - C^B a 

5. = 2Q,, (TC* = -Q<* 

6. V M Q„ = D^Qy - Q^Q U + \Q 2 gp V 

7. V a Q a = D a Q a + Q 2 , V a Q a = D a Q a + 2Q 2 

8. V a Q 2 = Q a Q 2 + 2Q»V a Q fl = <3 M V a Q M + Q„V Q Q" 

9. (V^V a - V a V*) Q„ = QpR^ a + Q»H afl = Q^R>\ 

We now define the auxiliary operator V a := V a — Q a ; thus V a commutes 
with g^ v . Note that V Q is not a derivative operator. For every torsion-free 
connection, the Bianchi identities take the form 

^ pR /ItJV Vjy-R ~\~ 0. 

We contract on a — o to obtain 

V ' pR^v — VuRpp + ^ a R a pu P = 0. 



14 



Transvection with g pp yields 

V^-V^+V^V = V^u-VuR-Va (g a "H^)+V a R^^ = 

2V^ V - V V R - V Q (g afl H^ u ) = 0. (A.26) 

On the other hand, 

R\ = g^R av = g pa (R (au) + H au ) =: R\ + g» a H au , 
hence equation (jA.26|) yields 

2V M i?^ - V„R + V a (g^H^) = 
and since V Q commutes with g atl we obtain 

V„ - ^Rj = -\v"H, u . (A.27) 

Denoting by G fl u the mixed tensor corresponding to the symmetric part of 
the Einstein tensor, i.e., 

the Bianchi identities (|A.27|) can be finally written as 

(V„ - Q M ) G\ = ~ (V - g^ 1 ) H^ u . (A.28) 

Note that in the case of integrable Weyl geometry, equation (jA.28|) reduces 
to V^G^V = Q^G 1 ^ which is usually referred in the literature as Bianchi 
identity (cf. [21). _ 

Relation with the Levi-Civita connection 

It is useful to express the Einstein tensor and the Bianchi identities in terms 
of quantities formed with the Levi-Civita connection, D. Starting with the 
definition of the Riemann tensor we arrive at 

o 

pa pa _j_ a ria a s-ya , /-ya sy[i sya 

n p^s — n -i- ^7^/35 — os^pnf t" — o^o^ 

I fa 1 I f/i \ /oa /a 1 f M \ /^a 

+ 1/nJ °/3<5 + 1/35/ °/n ~ °/37 ~~ 1 07 J 
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where the accent ° denotes a quantity formed with the Levi-Civita connec- 
tion. We contract on a, 7 and take account of = 2Q fl , g^C™ = 
—Q a , {a[i} = <9 M In y/—g to obtain for the symmetric part of the Ricci 
tensor, 

o xi 

R(ps) = Hps - DipQs) + -jQpQs - (D^QtM + Q 2 ) 

and upon a new contraction, 

R = R — 3D a Q a - hf. 

Therefore the Einstein tensor takes the form 

G (aP) = G Q(3 - D [a Q p) + ~Q a Qf, + g aP \D»Qp + ^Q 2 ^j . (A.29) 

We now turn to the Bianchi identities. We calculate the right-hand side 
(RHS) of (PT28|) 

(V„ - Q a ) G\ = D a G% + CysK - CLG a j_ - Q a G\. (A.30) 

Taking into account of (jA.29|) the first term in the RHS of (jA.30|) can be 
written as 

DgG a u = 

Dj a v -]pQ u -^D a D v Q a +D u D a Q a + l -Q v D a Q a + l - (Q a D a ) Q U +±D U Q 2 . 

We use the anticommutativity of the derivative operator D a to write the 
terms -\D a D u Q a + D u D a Q a as 

1 

—Q a R a u + -D a D u Q a . 
The second and the last terms in the RHS of (|A.30J) can be written as 

Ca f\ riot /-\ ria 

Q a G% - \ (Q a D a ) Q u - -D V Q 2 + l -Q v Q 2 + Q u D a Q a + ^Q U Q 2 . 



16 



The third term in the RHS of (jA.30|) can be written as 

-Ci v G% = - l - (8iQ v + %Q a - g au Q p ) G% = \q v R = 

1 ° 3 3 
-Q U R — -Q u D a Q a — -Q U Q ■ 

Putting all these together we obtain from equation (jA.30|) 

(V a - Q a ) G% = ~UQ V + l -D a D v Q a . 

On the other hand, 

(V Q — Q a ) H au = g a ^ (D^Huv — C? a H pv — C? u H ap ) — Q a H au = 
UQ V - D a D u Q a - Q a H au . 

Therefore the Bianchi identities (jA.28|) reduce to 

Q a H au = 0. (A.31) 
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